Spatial Extent of Random Laser Modes 
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We have experimentally studied the distribution of the spatial extent of modes, and the crossover 
from essentially single-mode to distinctly multi-mode behavior inside a porous gallium phosphide 
random laser. This system serves as a paragon for random lasers due to its exemplary high index 
contrast. In the multi-mode regime we observed mode competition. We have measured the dis- 
tribution of spectral mode spacings in our emission spectra, and found level repulsion that is well 
described by the Gaussian orthogonal ensemble of random-matrix theory. 

PACS numbers: 42.55.Zz, 42.25.Dd, 78.55. Cr 



Random lasers, in which gain is combined with mul- 
tiple scatteringof light [3, Q, are the subject of active 
research 0, 0, HI, El- The discovery of intriguingly nar- 
row spectral emission peaks (spikes) in certain random 
lasers by Cao et al. Ml has given an enormous boost to 
this field of research||, ©, M, [□], H H 0. In the 
past few years a debate on the origin of these spikes has 
started \U 
Cao et al. 



MM HI- On one hand 
as well as Sebbah and Vanneste and 
Apalkov et al. [17| . attribute the spikes to local cavities 
(local modes, LM) for light, which are formed by multiple 
scattering. On the other hand Mujumdar et al. [l4[, as 
well as Pinheiro and Sampaio [2l|, attribute the spikes 
to single spontaneous emission events that, by chance, 
follow very long light paths (open modes, OM) in the 
sample and hence pick up a very large gain. The LM 
and the OM model we cite represent divergent answers 
to the question what the spatial extent is of the modes 
responsible for the spikes. 

The spatial extent of the modes is a crucial factor for 
the fundamental behavior of a random laser. If the spa- 
tial extent of the modes is small, and the modes do not 
spatially overlap, the laser effectively consists of a col- 
lection of single-mode lasers. In contrast, spatially over- 
lapping modes inside the random laser lead to distinctly 
multi-mode behavior, such as mode competition (23J. 

In this Letter, we present the first systematic study of 
the spatial extent of the modes of a random laser, and 
show the crossover from essentially single-mode to multi- 
mode behavior in a new type of random laser: single 
crystal porous gallium phosphide filled with liquid 
dye solution. Among other conclusions, our results in- 
dicate that the LM model describes the physics of our 
random laser better than the OM model, but the LM 
model needs more sophistication. 

Our random laser consists of porous gallium phos- 
phide obtained by anodic etching (porosity, 45% air; 
thickness, > 46/im), of which the pores are filled with 
a 10 mmol/1 solution of rhodamine 640 perchlorate in 



methanol (pump absorption length, 22 /im; minimal gain 
length, 12 /im 25]). This random laser combines three 
characteristics which sets it apart from other random 
lasers. Firstly, the position of the scatterers is rigidly 
fixed. Secondly, the efficient laser dye rhodamine 640 
ensures high gain. Finally, the contrast of the refrac- 
tive indices in the sample is the highest ever reported 
for random lasers: n = 1.33 for methanol and 3.4 for 
gallium phosphide. The method of fabrication of the 
porous gallium phosphide is slightly modified from the 

HE?!- We used n-GaP with a 
obtained from Marketech. 



method described earlier 
dopant density of 2-10 17 /cm 2 



To make porous GaP without a non-porous top layer, 
the GaP is anodically etched in 0.5 M H2SO4. In the 
first step, the potential is set to 22.5 V until 5 C/cm 2 
has passed through the sample. In the next step, the po- 
tential is switched to 30 V for 10 minutes. In the third 
step, the potential is switched back to 22.5 V. The GaP 
is etched until 60 C/cm 2 has passed through the sam- 
ple. After finishing the etching process, the first layer 
can easily be removed, leaving porous GaP. 

The transport mean free path £ of our samples filled 
with pure methanol is measured with an enhanced- 
backscatter-cone experiment [28[ and is 0.6 ± 0.3 /im 
at A = 633 nm (k^t — 6.4 ± 2.8, with k$ the vacuum 
wavenumber, n e ff = 2.1 ± 0.4). The sample holder has 
a sapphire window with a thickness of 3 mm. The ran- 
dom laser sample is lightly pushed onto the window with 
a spring. A detailed calculation shows that the emit- 
ted laser-light is generated inside the pores of the gal- 
lium phosphide, and cannot be generated by the dye be- 
tween the window and the sample. The sample holder 
can be translated in one direction (parallel to the sur- 
face of the sample) with a total displacement of 12 mm 
by means of a stepper motor. The pump source is a 
tunable optical parametric oscillator (Coherent Infinity 
40-100/XPO). The pump pulse, at a wavelength of 567 
nm, has a duration of 3 ns and a rep rate of 50 Hz. The 
wavelength of the pump source is chosen to avoid dam- 
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FIG. 1: (Color) Two emission spectra of the random laser col- 
lected at the same spot on the sample. One spectrum (black) 
before, and one (red) after shifting the sample by 37 /im and 
back. Apparently, the frequencies of the modes are fully de- 
termined by the realization of the disorder. 
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FIG. 2: (Color) One of the 50 single-shot-emission spectrum 
of our random laser system collected at one position of the 
sample. The symbols near the peaks correspond to the pulse 
energy of those peaks in Fig. [3] 



age through absorption of light by gallium phosphide, 
and to obtain high enough gain of the dye. The pump 
light is spatially filtered, and focused with a microscope 
objective (numerical aperture, 0.55) onto the sample (fo- 
cus area, 3 ± 1 //m 2 ). The pump energy on the sample 
was at maximum 0.32 /xJ/pulse. Emission light of the 
random laser is collected through the microscope objec- 
tive, dispersed by a spectrometer (resolution, 1 cm -1 ) 
and detected by an electron-multiplier charged-coupled 
device (C9100-02, Hamamatsu) [35J. The responsivity of 
the system was calibrated with an Ophir power meter 
(PD300-3W) and a HeNe-laser. 

All measurements presented in this Letter are mea- 
sured at an input energy of twice the threshold input 
energy, except Fig. [3] This threshold is defined as the 
collapse of full width at half maximum (FWHM) of the 
emitted spectrum [2J. The FWHM of the spectrum far 
above threshold is 13 times narrower than the FWHM of 
the emitted spectrum below threshold. 

To discriminate between current models, and possibly 
stimulate new theoretical approaches, we want to inves- 
tigate whether mode frequencies are completely deter- 
mined by the realization of disorder inside the random 
laser. In contrast to the frequencies of the modes, we 
expect that in all cases the peak intensities will depend 
on the pump energy; as we have pulse-to-pulse fluctua- 
tions of our pump source, we expect a different overall 
scale factor for our peak intensities for every shot. We 
collected several spectra above threshold at a single posi- 
tion on the sample. In between measurements, the sam- 
ple was translated 37 /jm and back. Two of the resulting 
spectra are shown in Fig. [T] 36]. The spectral positions 
of the modes reproduce, proving that mode frequencies 
are completely determined by the realization of disorder, 
and not selected by spontaneous emission events. 

In addition to the expected overall scale factor the 
peak intensities above threshold differed individually 
from shot-to-shot, as is shown in Fig. [TJ To investi- 
gate the origin of this height-distribution variation, we 



looked at the emission spectrum of the random laser at 
a single position for different pump energies. We took 50 
single shot measurements, and used the measured pulse- 
to-pulse variation of the pump source to obtain a range 
of pump energies. Above threshold we observed a linear 
relation between the total energy of the emitted light and 
the pump energy. A typical emission spectrum is shown 
in Fig. [2] We determined the energies in the observed 
4 narrow peaks in the emission spectrum (integration of 
spectrum over 0.576 nm). In Fig. [3] we plotted these en- 
ergies versus the total pulse energy of the emitted light. 
When the total emitted pulse energy increases above a 
threshold of 50 pJ, we observed that a single mode starts 
to lase (star symbol). In this regime, the random laser 
can be described as a single-mode laser. We found a 
decrease of the spectral width of this single mode for 
increasing pulse energy in the mode, in correspondence 
with the well-known Schawlow-Townes behavior [29| . At 
higher emitted pulse energies (125 to 150 pJ) two more 
modes of the random laser start to lase, whereas the en- 
ergy in the first mode does not increase anymore. This 
leveling-off is a clear sign that the modes are compet- 
ing for the available energy, and thus are overlapping in 
space. This spatial mode competition causing the height- 
distribution variation is not related to the frequencies of 
the modes. 

To directly observe the spatial extent of the modes 
inside the random laser, we measured emission spectra 
while pumping at different positions of the sample. In 
contrast, using speckle correlation techniques it is pos- 
sible to determine the spatial extent of the secondary 
light source on the surface, originating from the scattered 
light (ljj . In our experiment the sample was horizontally 
translated with steps of 506 ± 4 nm, and for every po- 
sition of the sample 50 single-shot-emission spectra were 
collected. We displaced the sample 20 /*m in total, which 
is 10 times the diameter of the pump spot. In Fig. 0] 
zoom-in of collected spectra at 5 consecutive positions 
are shown 37[ ■ Figure H] is a visualization of a typical 
evolution of a peak in the emission spectrum as the po- 
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FIG. 3: (Color) Pulse energy in 4 modes [see for single-shot 
emission spectrum Fig. [2] versus the total emitted pulse en- 
ergy of the collected light for 50 single-shot spectra. The 
measured errors are smaller than the symbol size. The solid 
line is a guide to the eye for the star symbols. Clearly, one 
mode starts to lase (star symbol), before the other modes 
begin to lase at higher pump energies. 
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Spatial extent of a mode ((am) 

FIG. 5: Measured distribution of the spatial size of the modes 
determined with a displacement measurement. The dotted 
line is a marker for the diameter of the focus of the pump 
source on the sample. 
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FIG. 4: Enlargement of a part of the emission spectra where 
a peak is visible. Between consecutive spectra, the sample is 
translated by 506 nm. 



sition of the sample is changed. First a peak becomes 
visible. Additional displacement of the sample leads to 
an intensity going through a maximum. The increase 
and decrease of intensity corresponds to the overlap be- 
tween the gain region and the mode. The spatial extent 
of mode is defined by the difference in sample position 
where the mode rises above noise level and the position 
where the mode falls below the noise level (noise level is 
5%). In the specific case of Fig. [4] the spatial extent is 
2.5 ± 0.5 /im. 

In Fig. [5] the distribution of the spatial extent of all 
the modes is plotted, together with a marker of the di- 
ameter of the focus of the pump source on the sample. 
Almost 80% of the modes have a spatial extent smaller 



than the diameter of the spot size of the focus (2 zxm) . 
The other 20% of the modes are larger than the diam- 
eter of the spot size of the focus. For the LM (OM) 
model, we would expect all modes to have a spatial ex- 
tent much smaller (larger) than the focus diameter of the 
pump spot. Apparently, our measured distribution dis- 
agrees with the expectation of either of the two models. 
The selection criterium for a mode to lase is its dwell 
time [13]. Our observed limited spatial extent of lasing 
modes suggests that modes with long dwell times have 
small spatial extent. This connection has never been put 
forward before, and may be unique to three dimensional 
strongly scattering media [HI ■ 

We will now focus on the statistics of spectral spac- 
ings of the modes. For 11 spectra taken at widely spaced 
positions, we determined the spectral spacing between 
two adjacent peaks. These 11 spectra are statistically 
independent but equivalent. The only difference is the 
position on the sample where the data was collected. 
Therefore we combine the 41 spectral spacings of these 
11 spectra in one distribution. In Fig. [S]we plot this dis- 
tribution of the spectral mode spacing. Surprisingly, our 
observations show that the statistics are determined by 
level repulsion, in complete analogy to level repulsion in 



quantum mechanics [311 . [32j . A first experimental indi- 
cation of spectral level repulsion was shown by Cao et al. 

To describe level repulsion quantitatively we compare 
our observed statistics of the spectral mode spacings to 
the results of random-matrix theory. Although random 
matrix theory does not necessarily apply here, it does 
describe avoided crossings. The distribution of mode 
spacings according to the Gaussian orthogonal ensemble 
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FIG. 6: Measured distribution of the spectral mode spacing 
for 11 emission spectra. The solid line is a fit of Wigner's 
surmise fl}. The fit parameter is the mean mode spacing A, 
we find A = 16.2 cm" 1 . 

(GOE) is given by Wigner's surmise [II] 

P(x) ^Cxe-™ 2 '^, (1) 

where P is the distribution of the mode spacing x of the 
modes, A the mean mode spacing, and C a scaling fac- 
tor. In Fig. [6] we plot the fit of Eq. ^ (black line), with 
a mean GOE mode spacing A of 16.2 ± 0.9 cm" 1 . When 
considering all modes in a multiple scattering system the 
level statistics is known to be Poissonian (exponential) 
[34| . In this experiment we only observe those modes 
that lase and they show level repulsion. Apparently, the 
selection mechanism for random lasing causes the level 
statistics to change from Poissonian statistics to level re- 
pulsion. 

We have performed a systematic study of random las- 
ing in porous gallium phosphide. The frequencies of the 
modes in our sample are determined by the realization of 
disorder. We have observed a variation in the height dis- 
tribution, and attributed this variation to spatial mode 
overlap and gain competition. The distribution of the 
spectral mode spacings shows clear level repulsion and 
can be described with the Gaussian orthogonal ensemble. 
In addition the (small) spatial extent of lasing modes was 
observed. We hope that our observations stimulate more 
theoretical work. 
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